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We theoretically calculate the charge-supercurrent through a ballistic graphene junction where superconduc- 
tivity is induced via the proximity-effect. Both monolayer and bilayer graphene are considered, including the 
possibility of strain in the systems. We demonstrate that the supercurrent at the charge neutrality point can be 
tuned efficiently by means of mechanical strain. Remarkably, the supercurrent is enhanced or suppressed rela- 
tive to the non-strained case depending on the direction of this strain. We also calculate the Fano factor in the 
normal-state of the system and show how its behavior varies depending on the direction of strain. 
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I. INTRODUCTION 

Since the discovery of graphene 1 , many interesting phe- 
nomena have been predicted in the context of quantum trans- 
port in this material 3 ' 5-8 . It has been demonstrated in sev- 
eral theoretical and experimental works that the conductiv- 
ity of graphene monolayer junctions at zero doping level (the 
charge neutrality point aka Dirac point) displays a minimum 
value in the short junction regime 1 ' 2 5 ' 8 ' 9 . The physical rea- 
son for the minimum conductivity in this regime is the exis- 
tence of evanescent modes that can transport current over a 
finite length L of the system. A unique aspect of graphene is 
that these evanescent modes are predicted to generate pseudo- 
diffusive characteristics of the quantum transport properties 
even in ballistic graphene samples 5 . 

Bilayer graphene is a basic carbon structure and has at- 
tracted considerable attention recently. This system consists 
involved two coupled graphene layers with prominent charac- 
teristics such as pseudospin and variable chirality 15 . The chi- 
rality of the massless Dirac fermions in monolayer graphene 
is locked to the momentum direction and consequently lies in 
the plane of the sheet. In bilayer graphene, however, massive 
Dirac fermions with perpendicular chirality to the sheet plane 
may occur. Bilayer graphene features a band-structure sim- 
ilar to a semiconductor with parabolic bands with a tunable 
charge excitation gap 14 . 20 - 42 " 46 a gapless graphene bilayer is 
more stable compared to its gapped equivalence and thus oc- 
curs naturally. Bilayer graphene has also shown anomalous 
phenomena such as half-integer quantum Hall effect, mini- 
mum conductivity at zero energy and 2tt Berry phase. 

Very recently, the response of graphene to strain has been 
intensively examined. Several unusual effects have been un- 
veiled, including the generation of very high pseudo-magnetic 
fields of order 300 T 10 . An interesting question in this con- 
text relates to if strain imposed on graphene, whether it be 
mechanical or thermal in origin, can be used to control its 
transport properties". This issue is also motivated by the 
well-known fact that strain imposed on silicon-based devices 
can enhance their functionality. A mechanical deformation 
of a graphene sheet will invariably generate scattering centers 
which effectively influences the hopping amplitude, and thus 
suggests that the transport of Dirac fermions should respond 
to the presence of strain. It is known that strain may change 
the physical properties of nanotubes drastically 1415 ' 34-36 . The 



strain can be induced in graphene via several routes, including 
mechanically 37 ^ 11 . 

Based on this idea, we address in this paper a novel 
class of the Josephson graphene junctions with the capabil- 
ity to sustain tuneable charge-transport at the Dirac point 
by means of mechanically induced strain. To demonstrate 
this, we first solve the Bogoliubov de-Gennes equations 
both for a strained monolayer and bilayer graphene-based 
superconductor|normal|superconductor (S|N|S) junction with 
heavily doped S regions, as is experimentally relevant. We 
then derive explicit analytical expressions for the Andreev- 
levels and use these to obtain the phase-dependent supercur- 
rent I((f>) in the short-junction regime 12 . Both the critical 
current I c and the I c Rn product is investigated for a range 
of doping levels in the N region, including the charge neu- 
trality point. Above, Rm is the normal state resistance. Fi- 
nally, we calculate the Fano factor F in the normal (non- 
superconducting) state and show how this is influenced by the 
presence of mechanical strain in the system. 

To describe strained graphene, we adopt the model used in 
Ref. 13 for monolayer graphene and also consider a similar 
model for strained graphene bilayer junctions. Our findings 
show that for a zig-zag (Z)-strain (see Fig.l) the transmis- 
sion probability of evanescent modes near the charge neutral- 
ity point is suppressed, which influences both the conductiv- 
ity a and the supercurrent. However, when the strain is ap- 
plied along the armchair ( A) direction instead, the transmis- 
sion probability is enhanced and correspondingly influences 
charge-transport in the system. These results point towards 
new perspectives within tunable quantum transport by means 
of induced strain in a graphene mono- or bilayer. This find- 
ing might also be of relevance in the field of spintronics (val- 
leytronics), since the strain also affects the pseudo-spin of the 
chiral fermions in graphene ' 15 . 

The paper is organized as follows: In Sec. II, we present our 
theoretical approach and derive a general expression for the 
normal-state transition probabilities describing both strained 
monolayer and bilayer junctions. In this section, we also dis- 
cuss an experimental setup for detecting our predictions. In 
Sec. Ill, the Andreev subgap bound state energies for both 
strained monolayer and bilayer S|N|S Josephson junctions are 
obtained and the subsequent results are discussed. We finally 
conclude our findings in Sec. IV. 
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FIG. 1: The schematical setup of a strained graphene Josephson 
junction. The two superconducting electrode interfaces are located 
at x — 0, L. There is also a controllable gate voltage for tuning 
the concentration of carriers (not shown). Z and A stand for zig-zag 
and armchair strains while Oi represent displacement vectors of the 
three nearest neighbors C atoms in the strained graphene. 9 g is the 
strained angle between C-C junctions which can be either larger or 
smaller than the non-strained value 60°. 



II. THEORETICAL APPROACH AND MODEL 



(n.n.n.) in graphene can cause the Dirac cone to be tilted, this 
effect vanishes under the influence of strain of order ~ 20% 14 . 
In this regime, the generalized Weyl-Hamiltonian used here 
gives very good agreement with ab initio calculations. We 
note that it is also possible to model strained graphene by in- 
cluding a fictious gauge-potential A 28 . Motivated by the re- 
sults of Ref. 6 for bilayer graphene junctions, we here adopt 
the same model as the strained monolayer for the strained bi- 
layer whereas the interlayer hopping t± is left intact as the 
strain is applied in-plane. We emphasize that in this paper 
the same model for clean doped bilayer graphene region is 
adopted as Ref. 6 and consequently the trigonal warping ef- 
fects may be neglected. 

In both cases, we consider an S|N|S junction with s- 
wave superconducting electrodes. Previous works have been 
considered the Josephson effect in non-strained monolayer 
graphene 8 ' 22 ' 30 . We assume that the Fermi wavelengths sat- 
isfy A at ^> As, corresponding to heavily doped S regions. 
In this regime, we may ignore interface details and consider 
the following x-dependent superconducting order parameter 
as depicted in Fig.l (see Ref. 31); 



A{T)e l<t>l x < 
AO) = { < x < L 
A(T)e** r x > L 



(3) 



The two graphene systems considered in this paper 
(monolayer and bilayer) are modelled via the following 
Hamiltonians 6 ' 16-18 ' 20 (S and B stand for single- and bilayer 
graphene under strain, respectively): 



By substituting the Hamiltonians Eq.(l) and Eq.(2) into the 
Bogoliubov-de Gennes equation 
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(1) where /i is the chemical potential, we find the following 
energy -momentum dispersion relations: 
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where ir=v x p x +iv y p y and v x . y are the Fermion velocities 
in the x, ^-directions while a XyV are Pauli matrices. Here, 
U represents an external gate potential. The ± signs refer 
to the K and K' valleys of monolayer graphene. We here 
consider an A2B1 stack for the bilayer graphene sheet and 
use a tight-binding model in which the massive chiral Dirac 
frmions are governed by Eq. (2). The resulting Hamilto- 
nian becomes similar to a gapless semiconductor Hamilto- 
nian with parabolic electron and hole bands touching when 
t— >oo. In order to model applied strain to the system, we 
adopt the model of Ref. 13 and expand the tight-binding 
model band structure with arbitrary hopping energies £1,2,3 

i.e. e — ± I Yli=i ti e ~ l k °^ I around the Dirac point, = 
(cos- 1 (-l/2ri)/V~3a Xl 0) H ' 24 . As shown in Fig. 1, at are 
displacement vectors between three nearest C atoms (see Ref. 
23). We assume tx,2 = t\- and t% = t in our calculations and 
set 77 as the ratio t\- / 1. This assumption generates asymmet- 
ric Fermion velocities along the different directions. These 
velocities are given by v x — 2t\-a x sm(cos~ 1 (—l/2r)))/h 
and v y = 3ta y /2h 14 ' 15,24 . The next-nearest neighbor hopping 
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u = v x cos + v y sin 



Above, 9 is the angle of incidence of the particles. The eigen- 
functions of Eq.(4) for the strained monolayer and bilayer sys- 
tems within the normal region are found to be: 
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where s and e represent 1x4- and 1 x 2-spinors with only 
zeroes as entries, while e and h stand for electron and hole 
particles. The ± sign refers to right and left going particles. 

TV 
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The coefficients a N ± and a\f+ are defined as follows 
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For the bilayer system, we consider the hopping value be- 
tween the graphene layers to be smaller than the doping level 
of superconducting regions, still much larger than the super- 
conducting gap that is /i s t ef,s . The former as- 
sumption assures ignoring the contact details in the S|N in- 
terfaces while the latter not only helps to simplify theoretical 
approach but also warranties realistic approximations in our 
analytical calculations. Within the normal region in which su- 
perconducting order parameter A=0, the Eq. (4) leads to un- 
coupled equations. In this paper, we focus on the low-energy 
regime. One then finds the following parabolic dispersion re- 
lation for electrons and holes in the normal bilayer region: 
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Due to translational symmetry in the transverse direction, k y 
and e are both conserved upon reflections at the interfaces 
located at x = 0, L. Accordingly, the dispersion relations 
and the following equation assure both energy and momen- 
tum conservation of particles in the y-direction upon Andreev 
electron-hole conversion, 
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The total spinors in the three regions thus read: 
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In superconductor spinors, we define the following relation 
for superconducting coherent factors (see Appendix): 
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the definition helps to simplifying our notation. The spinors 
in the superconducting regions carry S superscript, whereas r 
and I stand for right and left superconducting regions. The su- 
perconducting phases in each region are assumed to be <p r and 
4>i, while th and t e are the scattering amplitudes of hole- and 
electron-like quasiparticles. Matching the total wavefunctions 
at each of the interfaces, i.e. 
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leads a quantization relation between superconducting phase 
difference <f> = </> r — <pi and the quasiparticle excitation en- 
ergy e. The boundary conditions lead to a 8 x 8 matrix M for 
transmission and reflection coefficients which is presented in 
the Appendix for monolayer case 32 . det(M) = generates a 
non-trivial relation between e and d> as follow; 



Fi + F 2 sin 2(3 + F 3 cos 2(3 = 0, 
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We here employ the most relevant experimentally approxima- 
tion i.e. the "short-junction" regime in which AL/hv <C 1. 
Within this regime, F 1 , F2 and F 3 are reduced to the follow- 
ing expressions; 



Fi = -« v _ + <+r cos + sin"(£;F) 

x((af + ) 2 + af-(-(af + ) 2 af-+« e V -)-l) 

F 2 = (16) 

F 3 - sin 2 (fcL)(af-af + + l) 2 + (a£ + af + ) 2 cos 2 (/fcL). 

In this case, using the definition of (3 and Eq.(15) the single 
Andreev bound state is obtained vs F 1 and F 3 
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and then 



£„(<£) = A J 1-T n sin 2 0/2 



in which r n is transmission probability for the normal 
graphene region between two strongly doped electrodes (ei- 
ther superconductor or normal). After some calculations we 
reach at an expression for r n valid for both strained monolayer 
and bilayer systems i.e. 



n (af + + af_) 2 cos 2 (fc„F) + (af + <_ + l) 2 sm 2 (fc„L) ' 

(17) 

We utilize the general r„ for investigating the transport prop- 
erties of the strained monolayer and bilayer junctions in the 
next section. 

The contribution of the Andreev bound-state spectrum to 
the supercurrent is given by 8 ' 12 : 



I(4>) = —T- 22 T n Sin (j}/£n{4>)- 



(18) 



n=0 



For a wide graphene junction, W^L, the boundary condi- 
tions (zig-zag and armchair) at the edges y—W/2 and — W/2 
are irrelevant and we here assume smooth boundaries for the 
two edges. In this regime, we may replace the summation 
over quantized modes with an integration: *}2 n -^W '/ir J dq n . 
Now we proceed in the following sections to study the super- 
current using Eqs. (8), (9) and (18) for monolayer and bilayer 
Josephson junctions in particular at the charge neutrality point 
i.e. fx 0. 
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FIG. 2: The critical current 7 C (left panel) and its product with the 
normal-state resistance I c Rn (right panel) as a function of \iLjhv 
for a monolayer system. The solid line pertains to a non-strained 
junction. For Z-tension tt- = 0.56to, t = l.lio while for yl-strain 
t\- — 0.95to, t — 0.5to- The arrow indicates how the Josephson 
current may be enhanced by means of the applied direction of ten- 
sion to the system. The inset panel shows critical supercurrent as a 
function of strain. 



FIG. 3: The critical current I c (left panel) and its product with the 
normal-state resistance I c Rn (right panel) as a function of fiL/hv 
for a strongly coupled bilayer system. The solid line pertains to a 
non-strained junction. Here, the same values as for the monolayer 
system have been used for the tensions, i.e. in the case of Z-tension 
t\- — 0.56to, t = l.lto while for A-tension t\- = 0.95to, t = 
0.5t fors = 0.2. 



III. SUPERCURRENT, FANO FACTOR AND ANDREEV 
BOUND STATES IN STRAINED GRAPHENE 
MONOLAYER/BILAYER S N S JUNCTION 

By inserting Eq.(8) for the monolayer system into Eq.(17), 
the following expressions are obtained: 

s = ["* cosflf] 2 

" [v x cos 9 s n cos(fc5 J L)]2 + [ v s sin(fc5L)]2 

k s - ^ 0-*m( h2vlql V 

In the case of a bilayer-system where Eq.(9) holds, the corre- 
sponding normal-state transmission probability takes the from 

T B = M [2«xCOSflg]7^ 

» M [2^cos^cos( fc gL)]^ + tx[Wtx + 1} sin (fcB J L)]2 

The behavior of / c and i c i?jv for strained monolayer 
graphene-based Josephson junction vs /iL/hv is shown in 
Fig.2. The critical current against strain is shown in the in- 
set panel (see Ref.25). For ^4-strain, we assume t\- = 0.95io, 
t = 0.5t (to — 2.7 eV for non-strained graphene) and 
for Z-strain t\- = 0.56£o, t = l.lto which follows when 
s = 0.2°. As seen, the critical current becomes nearly zero at 
the charge neutrality point for Z-strain, whereas the current is 
enhanced compared to the non-strained case when the strain is 
applied along the A direction. In effect, the Z-strain induces 



a very small phase-dependent contribution to the Andreev- 
bound states and the supercurrent vanishes. This suggests a 
remarkable fact: the supercurrent can be efficiently tuned by 
means of both the magnitude and the direction of the strain 
imposed on the system, even at the Dirac point. We have also 
considered the same model of strain for a bilayer graphene 
S|N|S junction and plot I c and I c Rn as a function of pL/hv 
in Fig. 3. In this case, it is seen that the critical current tends 
toward zero as one approaches the Dirac point because of 
the assumption t± S> e, A which influences the evanescent 
modes 26 . The obtained normal-state transmission probability 
t® is proportional to p and hence tends toward zero as one 
approaches the charge neutrality point. In order to understand 
these results, one has to consider two facts: i) all transport 
modes n in the system become evanescent (k n = iq n ) at the 
Dirac point and ii) have a transmission probability through the 
junction given by r„. In the Z- and yl-strain cases, t„ decays, 
respectively, faster and slower than the non-strained graphene 
as a function of q n . In turn, this dictates the magnitude of the 
contribution of transverse modes to the electron transmission 
and thus to the discrete Andreev bound state spectrum for the 
A- and Z-strains with respect to the non-strained system. We 
have also calculated the Fano factor (the ratio of noise power 
and mean current) via the normal transmission probability r n , 

defined as F = 2~^o° T n (l — 7"n)/X)o° r ™ 5 - ^ ne resu l ts f° r 
both mono- and bilayer graphene with and without strain are 
shown in Fig. 4. In the non-strained case, we reproduce pre- 
vious results for monolayer 5 and bilayer 6 junctions where a 
weakly doped middle region is sandwiched between two heav- 
ily doped regions 7 ' 27 . The scenario with strain has not been 
considered up to now, and inspection of Fig. 4 reveals that the 
strain influences how F evolves with the doping-level in the 
middle region. More specifically, in the Z-strain case the con- 
tribution of the transversal modes is suppressed and therefore 
the F goes towards saturation faster than the A- and non- 
strained regimes as the doping-degree p is increased. The sys- 
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FIG. 4: The Fano factor F for both monolayer and bilayer systems 
with a weakly doped region sandwiched between two heavily doped 
sides (the normal-state of the mentioned S|N|S junction) as a func- 
tion of nL/hv. In the scenarios with strain, the same parameters as 
in Fig. 2 and Fig.3 have been used. 



tem under tension, however, sustains still the universal value 
of F = 1/3 at the Dirac point just as the non-strained mono- 
layer system or diffusive normal metal 5 6 . We note that the in- 
fluence of trigonal warping may be neglected in the monolayer 
case when the impurity-potential is weak (ballistic regime) 29 . 
For the bilayer case, the trigonal warping becomes influen- 
tial in the low-energy regime |e| < 0.57i(t>3/i;) 2 where a 
relevant estimate for the parameters is 71 = 0.39 eV and 
vs/v = 0.1 -C l 17 - This yields |e| < 2 meV. However, the 
influence of strain in the considered bilayer model in this pa- 
per becomes most evident at higher doping levels as seen from 
Fig.3 where the trigonal warping effects can be neglected. 



IV. CONCLUSION 



mal transition probability valid for both strained monolayer 
and bilayer graphene systems. We have demonstrated that the 
direction of the applied strain to the system near the charge 
neutrality point can be used to efficiently tune the magnitude 
of the supercurrent in such a system. In addition, we have con- 
sidered the Fano factor F in the normal-state of this junction 
and how it is influenced by strain in the system. In this case, 
we also find that the direction of the strain is influential with 
respect to how F depends on the doping-level of the graphene 
sheet. We believe that these results point towards new per- 
spectives within tunable quantum transport by means of me- 
chanically induced strain. Interesting phenomena may be ex- 
pected to arise out of the coexistence of proximity induced 
ferromagnetism and superconductivity in a strained graphene 
junctions 33 . 

Appendix A: Andreev subgap states 



In this appendix, we present more details of our analytical 
approach used to find the general normal transition probability 
in Sec. II. We here also examine our analytical expressions 
for non-strained monolayer case where s is equal to zero. In 
the graphene monolayer superconducting regions, the right- 
and left-going quasiparticles are described via the following 
spinors 
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Similar spinors are obtained when starting with the Hamil- 
tonian Eq. (2) for the strained bilayer case, although they 
become 1x8 arrays. We focus our attention on the strained 
monolayer Josephson junctions in this appendix. Matching 
the wave functions of the superconducting and normal seg- 
ments at the two interfaces generates the following matrix for 
reflection and transmission coefficients. 



In conclusion, we have proposed a novel class of ballistic 
graphene monolayer/bilayer-based Josephson junctions with 
mechanical strain. We have derived a general analytical nor- 
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To determine the relation between the subgap energy of the use the determinant of M as mentioned in the Sec. II. Previ- 
quasiparticles and the superconducting phase difference, we ously in this paper, we have considered a heavily doped su- 
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perconducting regions which dictates normal trajectories of conducting region i.e. fi s > e, A and then Of ~ Of — 7 7^ 0. 
the quasiparticles relative the interfaces inside these regions. In this regime, we find Fi, F2 an d F 3 factors as follow; 
In this appendix, we now allow for a moderately doped super- 



F 1 = sin(fc^L) sin(fc^ L)(sin(7) — sin(#))(sin(7) + sin(^)) — cos 2 (7) cos(#) cos ( 0,4 ) cos((/>) 
F 2 = cos(fcf L) cos( 7 ) cos(0) sin(/c^ J L)(sin(7) sm(0 A ) + 1) 
— cos(fc^ L) cos(7) cos{9a) sm{k^ L)(sin(7) sin(0) — 1) 

F 3 = cos(fcf L) cos(k^L) cos(0) cos(0^i) cos (7) + (sin(7) sin(0) — l)(sin(7) s'm(0A) + 1) 



We denote e = 6 and 9/ t = 9 a- If we apply the short junc- superconducting regions i.e. 7 — > 0, we recover the results of 
tion approximation to the factors and assume heavily doped Ref. 8. 
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